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We examine a scheme for the optical creation of a superfluid vortex in a trapped Bose-Einstein conden- 
sate (BEC), using the stimulated Raman adiabatic passage (STIRAP) technique. By exposing an oblate, axis- 
symmetric condensate to two co-propagating laser pulses, one can transfer external angular momentum from 
the light field to the matter wave, if one of the beams is the fundamental Gaussian mode and the other is a 
Gauss-Laguerre mode of angular momentum lh. We demonstrate the complete transfer efficiency by numeri- 
cal integration of the multi-component Gross-Pitaevskii equation and explain the results with an intuitive and 
accurate approximation within the Thomas-Fermi limit. In addition, we discuss residual excitations (breathing 
modes) which occur in the two-dimensional regime and present the Bogoliubov excitation spectrum. 
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I. INTRODUCTION 

Ultra-cold atomic gases have provided us with novel phys- 
ical systems that exhibit all genuine many-body phenomena 
known from traditional condensed matter physics and, still, 
admit all the superior coherent control tools used in quan- 
tum optics. After the experimental realization of BEC itself, 
tremendous efforts were focused on the creation of topological 
and solitary excitations of condensates (for a current review 
see Ref. [l]]). Especially alternative methods for the creation 
of vortices have stirred the minds, as the traditional "rotating- 
the-squeezed-bucket" procedure was not successful, initially. 
Thus, the first fruitful proposal to create a vortex involved a 
rapidly rotating Gaussian laser beam entangling the external 
motion with internal state Rabi-oscillations |2, 3]. Later, con- 
densates were stirred up mechanically EttH. lal an d evapora- 
tive spin-up techniques created vorticity (7j and recently giant 
vortices 1 8] could be created. Due to larger asymmetries in the 
trapping potentials that can be achieved nowadays, vortices 
are now predominantly created with the stirring method and 
fascinating Abrikosov lattices containing up to 300 vortices 
have been made Alternativel y, p rospects for creat- 

ing vortices by optical phase imprinting II ill were investigated 
(in analogy to the successful soliton experiment |12jn and ap- 
plying magnetic interactions were considered 11311 . More- 
over, there have been ideas to create vorticity by sweeping 
a laser beam on a spiraling trajectory across the trap, induc- 
ing a Landau-Zener transition between the irrotational and the 
rotational state 11411 . 

On the other hand, the transfer of angular momentum from 
an optical field to a macroscopic rigid body or an atomic par- 
ticle has also a long standing tradition in quantum mechan- 
ics. The first proof that circularly polarized light carries an- 
gular momentum dates back to Beth's original experiment of 
1936 1 15]. More recently, due to the availability of Gauss- 
Laguerre laser beams with well defined external angular mo- 
mentum flrjl [l7ll it is possible to use them as optical tweez- 
ers and twisters 1181 U9H2OII . Surprisingly, even the transfer 
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of angular momentum to ultra-sonic waves in fluids can be 
achieved that way EUl . In the context of a BEC, using the 
angular momentum of light to create a doubly charged vortex 
has been proposed ll22ll . 7r-pulses in Raman type transitions 
were examined l23ll and an adiabatic passage to a vortex state 
was investigated by changing the two-photon detuning of an 
effective two-level system 12411 . 

In this paper, we will examine the transfer of external an- 
gular momentum of light to the matter wave with the help of 
a stimulated Raman adiabatic passage (STIRAP) l25ll . The 
basic effect relies on a quantum mechanical interference be- 
tween two ground states and gives rise to a multitude of phys- 
ical phenomena, e.g., dark resonances in optical spectroscopy 
pal, velocity selective coherent population trapping (VSCPT) 
12711 . a drastic modification of the optical index of refraction 
(EIT) of normal and BE condensed systems I29L 13011 . as 
well as constructive procedures to prepare l3lll and readout 
quantum states of atomic beams and optical cavities ll32ll . 

This article is organized as follows: In the Sec. II we will 
develop a scheme for creating a vortex in a BEC using the STI- 
RAP method. In analogy to single-particle physics, it is pos- 
sible to derive the relevant three-component Gross-Pitaevskii 
equation. In Sec. Ill we will present the results of numerical 
calculations that are in good agreement with a simple ana- 
lytical approximation within the Thomas-Fermi limit. In ad- 
dition, we will discuss the physics of the remaining residual 
excitations in terms of the "breathing modes" of a quasi two- 
dimensional system. In particular, we will calculate the radial 
Bogoliubov excitation spectrum of a condensate in the vortex 
state with angular momentum lh. Finally, we will summarize 
our results and conclude in Sec. IV. 



II. STIRAP IN A BOSE-EINSTEIN CONDENSED GAS 

The STIRAP method is now applied to a trapped BEC of 
three-level atoms in a A-type configuration shown in Fig. [2 
The two internal electronic ground states, e.g. the hyperfine 
levels of an alkali atom I29H33H34I1 . are denoted by \b) and |c), 
and by absorbing an optical photon, one reaches the excited 
state I a), respectively. The condensate is confined spatially by 
an oblate, axis-symmetric harmonic potential, which reads in 
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cylindrical coordinates (r, (j>, z) 



1 



y(r,z) = -M( W V+^V), 



(1) 



where M denotes the single-particle mass. By choosing the 
radial trapping frequency much less than the longitudinal fre- 
quency, i. e., uj <C u> z < one can confine the motion effectively 
to the radial component. 
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FIG. 1: A A-system exposed to two co-propagating laser pulses 
in two-photon resonance. The atomic transition frequencies are 
Doppler- shifted with respect to the rest frame i.e., uj a = cu a + 
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Laguerre laser beam (td) carries 1ft of angular momentum. The in- 
dividual detunings from the excited state are A p = uo a — (jj b — cu p , 

Ad = LOa — U>c — Ud- 



Now, we will expose the dilute atomic gas to two co- 
propagating traveling "pump" and "dump" laser pulses, 



E p (z,t) = \ p e p (t)e 



-i(ujpt — k p z) 



(2) 



E d {r,ct>,z,t) = \ d e d {r, 4>,t) e^"**-^ +c.c, (3) 

where A p and denote the corresponding polarization vec- 
tors. The slowly varying laser beam envelopes e p , Ed have a 
non-trivial temporal and spatial structure 



s p (r,t) w £p (t) = e p 



- P p -{t-r) 2 /2d 2 



£d(r,<j>,t) 



e d re' r e 



i<t> p ~t 2 /2d 2 



(4) 
(5) 



For the pump pulse, we choose the fundamental Gauss- 
Laguerre (GL) laser mode fl6ll with a spatial width larger 
than the BEC size and a temporal Gaussian turn-on shape 
with the width d. Numerical simulations prove that in this 
limit the spatial Gaussian envelope can be disregarded alto- 
gether in Eqs. 1415 1 . This pulse reaches its maximum inten- 
sity at some time t > 0. In order to transfer orbital angular 
momentum from the light beam to the matter-wave, we pick 
the first excited GL mode that carries external angular mo- 
mentum, the so called "doughnut-mode", for the dump beam 
fl7l f]~8L I22L l23ll . While spatial extension and temporal du- 
ration d can be set equal in both pulses, it is crucial that the 
dump beam reaches its maximum intensity at t = 0, first. 



This "counter-intuitive" pulse sequence is the key of the STI- 
RAP procedure and achieves an efficient adiabatic passage for 
linear, dissipative quantum systems I25I l27t l32ll . A full pop- 
ulation transfer can be reached if the field amplitudes satisfy 
the conditions 
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(6) 
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which is guaranteed by the pulse sequence given in Eqs. 

m 

Far below the transition temperature T <C Tbec, one can 
describe a multi-component BEC effectively within a simple 
mean-field picture 1351 l3ol . Thus, we introduce a three com- 
ponent state vector ^(r, <), that represents the components of 
the macroscopic atomic matter- wave. For the time evolution 
of this multi-level state vector, one can derive a generalized 
Gross-Pitaevskii (GP) mean-field equation l37il 



*(r,i) = (* a (r,i),* b (r,t),* c (r,t)) T , 
ihd t ^(r,t) = H(t)V{r,t). 

\I/(r, t) is normalized to the total particle number 



(8) 
(9) 



N = J d 3 r(|* Q (r,t)| 2 + |* b (r,i)| 2 + |* c (r,i)| 2 ), (10) 

in the BEC. Due to the unitarity of Eq. l|9}, this particle num- 
ber is conserved at all times. However, as we use explic- 
itly time dependent laser fields, the energy of the system can 
change (see Sec. Ill B). 

The internal structure of the Hamiltonian is quite easy 
to understand, as it follows straight from the single-particle 
physics that rules the dynamics of the dilute gas interacting 
with light. In Fig.^ we have depicted the optical dipole tran- 
sition scheme for a A-type atom. Within the standard rotating- 
wave approximation of quantum optics l39ll . one finds for the 
internal state Hamiltonian 



H(t)/h = ( n;(t) h b + s 
,n;(r,<M) h c -s 



(ii) 



The Rabi frequencies il p (t) = e p (t)d ba /h and analogously 
Od(r, c/>, t) = Ed{r, <fi, t)d ca /H, measure how well the photon 
field couples to the electronic transition and are proportional 
to the atomic dipole moments d ba , d ca . The remaining pa- 
rameters are the Raman detuning A = (A p + A^)/2 and the 
two-photon detuning S = (Ad — A p )/2, which refer to the 
individual detunings A p , A^, of laser the frequency and the 
Doppler-shifted electronic transition frequency. In order to 
achieve the optimal STIRAP performance, we will assume a 
two-photon resonance condition S = later, and pick a non- 
vanishing detuning A in order to avoid detrimental sponta- 
neous emissions, which would disrupt the coherent evolution. 
In the preceeding derivation of Eq. il It . we have also tacitly 
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adopted co-moving and co-rotating reference frames such that 
*(r,i) = (e l[(i '"- A) *-^ 2l * a (r,i),e l(i '^ ,5) ** b (r,t), 

el [(i, c +5)t-(fe p -/c d ) z ]^^ rji ^ T ; (12) 

in order to strip off the trivial plane wave character of the 
beams, see Eqs. i l2!3t . To keep the notation simple, we will 
drop the bar over the state vector in the following discussion. 

For the external motion of the multi-component gas, we 
want to assume that all components are confined by the same 
harmonic potential, cf. Eq. Q. This is not a stringent require- 
ment for the excited state component ty a , as particles will only 
rarely occupy this level and move through it quickly. Due to 
this low occupancy, it is also not necessary to consider any 
mean-field shifts arising from self-interaction, or the motion 
through the remaining components. Consequently, we can 
simply use the bare trap Hamiltonian with a Doppler shift 
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For the ground state components, we will only consider mean- 
field shifted energy contributions that arise from elastic colli- 
sions and denote the inter- and intra-species scattering lengths 
by {abb, o-bc, o-cb, 0-cc] ■ Hence, these components read 



hhb = wi + y(r ' z) + ^w {ahb |v " fo|2 + abc |v " c|2) ' 



Alth 2 , , _ ,2 I 1 2 \ 

-^(acci*cr+ocii*6r). 



(14) 



(15) 



The optical absorption-emission cycle imparts angular, as 
well as linear momentum onto the final state matter- wave f c . 
However, linear momentum cannot be conserved in a trapped 
system, and it would lead to a sloshing motion along the z- 
direction. This can be suppressed either, by choosing equal 
laser frequencies and photon momenta, or by squeezing the 
trapping potential into a very oblate configuration such that 
(3 = u> z /u> ^> 1, This effectively "freezes" the longitudinal 
motion due to an energy selection argument. 

The later situation leads to a more stable configuration and 
is favorable. Thus, we are able to approximate the state vector 

(16) 

by factorizing it into a radial part, which is normalized to the 
number of particles in the BEC, i.e. 

/>oo 

N = 2tt drr (|^ a | 2 + |^ & | 2 + |^ c | 2 ) , (17) 
Jo 

and the ground state 

Vo(«, t) = (/3/^) 1 / 4 e-^/ 2 e-^ 2 / 2 (18) 



of the one-dimensional harmonic oscillator in z-direction, 
which is normalized to one, i.e. 



dz|<Po| =1- 



(19) 



Please note that the third component now carries the mechan- 
ical angular momentum of 1ft per particle. After project- 
ing the state vector along the z-direction, we finally arrive at 
an effective Hamiltonian in the radial direction. Using the 
natural units of the transverse harmonic oscillator (time unit 
Tho = 2-k/lj, length scale oho = \J h/Muf), one finds 



H(t) 



id t iP(r,t)=H(t)ip, 

a n p (t) 



(20) 



'ft(°)+A Op(t) rn d (t) 

n*(t) h^+S + K n{r) 
rfl* d (t) ftW - S + Kn(r)j 



(21) 



where k = y/SnP a and all scattering lengths are equal to 
a. This assumption holds well for 87 Rb (see below). The 
particle density is denoted by n(r) = | - 0f)( r )| 2 + IV'cMI 2 , an d 
hS m ' represents a two-dimensional radial harmonic oscillator 
Hamiltonian 



r 



(22) 



r r* 

in an angular momentum manifold with m 0. 

III. RESULTS AND DISCUSSION 

In the following, we will be more specific and choose the 
typical scattering parameters for 87 Rb ll33l l34ll . To simplify 
matters, we will also assume that the self-scattering and cross- 
component scattering lengths are equal: abb = a cc = a& c = 
a = 110 Bohr radii (ao). This turns out to be a very robust 
approximation and we comment on this later. By implement- 
ing a numerical algorithm for solving the three-component 
Gross-Pitaevskii equation Eq. i20i , we can show that a quan- 
tized vortex (m = 1) is building up. For the trap frequen- 
cies we choose w = 2ir ■ 10 s" 1 and lo z — 2ir ■ 1000 s -1 , 
respectively. With an atomic mass of M = 86.91 amu, we 
find a harmonic oscillator size ano = v/ h/Mu> = 3.41 /im. 
In this spatial unit, we get for the scaled coupling constant 
k = ^/8Tr/3a = 0.0855. We will consider an atomic ensem- 
ble with N = 10000 particles. From a simple Thomas-Fermi 
approximation (see Eq. ( 1291 ) one finds estimates for the con- 
densate radius Rjp = = 5.7 with a chemical poten- 
tial fijp = 16.5. The laser parameters are chosen as follows: 
r = 0.3, d = 0.15, ft p = 200, fl d = 200, A = 30, and S = 
(two-photon resonance). In practice, A will be of the order 
of at least 10 5 , but for the simulations much lower values are 
sufficient. 

In Fig. 2 the transfer efficiency is demonstrated. Initially 
the whole system is prepared in the ground state \b). While 
the excited state \a) is not populated significantly at any in- 
stant, the population of state \c) is rising up. An almost 100% 
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transfer to the vortex state is possible and the total number of 
particles in the BEC is conserved as required by unitarity. The 
time evolution of the vortex state density is shown in Fig. 3. 
At some intermediate time (t = 0) high-frequency excitations 
appear. After the transfer is completed (t = 1), the shape of 
a vortex with angular momentum lh is matched almost per- 
fectly with the stationary single charged vortex solution. Still, 
some minor excitations, i. e., the "breathing mode" remains, 
which will be discussed in Sec. Ill B. In Fig. 4, we demon- 
strate the time evolution of the density of the irrotational state 
| b). Initially, the atoms of the condensate are prepared in the 
ground state (t = —1). This population is decreasing during 
the STIRAP process and vanishes afterwards. 



N(t) 
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FIG. 2: Time evolution of the individual level populations vs. time 
measured in natural units of Tko- The dotted line depicts the pop- 
ulation N a (t) of the excited state, the dashed-dotted and the solid 
lines represent Nt(t) and N c (t), respectively, and the dashed line 
corresponds to the total number of particles in the system, which is 
conserved, i.e. N(t) = N a (t) + N b (t) + N c (t). 



ficiency. It only leads to a further increase of the "breathing" 
of the vortex. As an example, we depict the results of the 
simulations for a vanishing cross-component scattering length 
(abb, a>bc, a cc ) = 104 ao (1.03, 0, 0.97) in Figs. 5 and 8. Once 
more, these simulations confirm the remarkable robustness of 
the STIRAP scheme. 
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FIG. 3: Time evolution of the vortex state density \ip c (r, t)\ 2 vs. ra- 
dius r measured in units of ciho. The solid lines depict the vortex 
state at an intermediate instant (t = 0) and after the STIRAP pro- 
cess is completed (t — 1). In general, the instantaneous Thomas- 
Fermi approximation (dashed line) compares well with the exact re- 
sults apart from minor excitations. 



A. The two-component Thomas-Fermi approximation 



While we have considered equal scattering lengths occur- 
ring accidentally in 87 Rb, this is not exactly valid and also 
not the generic situation applicable to other condensed ele- 
ments. It could be imagined that the imbalance of mean-field 
shifts deteriorates the transfer efficiency due to a violation of 
the bare two-photon resonance condition. Still, one could ap- 
ply a time dependent detuning that rectifies this effect in order 
to achieve an optimal adiabatic passage. Such a strategy has 
been proposed in Refs. I22LI24I1 . 

However, we have studied numerically the situation of un- 
equal scattering lengths without modifying the detunings and 
found no averse effects. First, we have examined the effi- 
ciency of the transfer procedure using the real scattering data 
of the JILA experiment I33I l34ll for the trapped hyperfine 
states \b) = \F = l,m F = -1) and |c) = |2, 1). Using the 
scattering parameters (abb, o>bc, o-cc) = 104 ao (1.03, 1, 0.97), 
there is an almost complete transfer to the rotational state pos- 
sible, yielding essentially identical curves as shown in Figs. 2- 
4. However, there are slightly higher excitations than in the 
idealized case. Second, we find that even a drastic varia- 
tion of the scattering lengths does not spoil the transfer ef- 



The Thomas-Fermi (TF) limit is an extremely useful ap- 
proximation for the ground state of an interacting BEC. Where 
applicable, it simply disregards the kinetic energy contribution 
to the total system energy. This is also possible for an excited 
vortex state. We will now consider this TF-approximation for 
the two ground-state components, after adiabatically eliminat- 
ing the ip a component, 



fi P (t) 



i> a (r, t) ps --^-^u 



0c 



(23) 



in order to obtain a simple estimate. During the adiabatic pas- 
sage, the irrotational state ipb shall be transformed coherently 
into the vortex state ip c . Thus, the relevant Hamiltonian in the 
reduced two-state manifold reads 



H TF (t) = — + nn(r)- 



n* d {r,t)n p (t) 

A 



Q, d (rd)Q,' p (t) 
4 



\n d (r,t 



1 

2r- 



(24) 

Once more, we want to assume an adiabatic following 
condition to obtain a stationary solution of ip(r, t) = 

cxp -i f ^drHTF^) (ifj b (r),ip c (r)) T - By solving the 
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simple ensuing eigenvalue problem of Eq. i24l . we find for 
the chemical potential in the vortex branch that 



(J-TF 



+ kutf + S/j,, 



2 r z 



(25) 
(26) 



|Oi| 4 -2|Oi| 2 



To keep the notation simple, we have dropped all the spatial 
and temporal arguments. Vice versa, one has to determine 
the chemical potential [It f such that N = 2ir dr r utf 
is satisfied at each instant of the adiabatic passage. For the 
vortex component of the corresponding eigenvector, one finds 



*JnTF 



COS( 



tan( 



smt 



A(5F2 - W-Sfi) 



(27) 
(28) 



It turns out that the chemical potential is almost constant dur- 
ing the transfer process. Hence, we can choose the value of 
Utf for the TF-solution of the vortex state, i.e. 

r 2 1 

Utf = HTFi^p.d = 0) « — + + nn TF = const. (29) 
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FIG. 4: Time evolution of the irrotational state density \Tpb(r)\ 2 vs. 
radius r measured in units of <2ho- The solid lines depict the posi- 
tion density at the initial (t = —1) and at some intermediate time 
(t = 0). The instantaneous Thomas-Fermi approximation (dashed 
lines) matches the exact solution after the transfer of particles starts 
effectively (t > 0). 



the vortex state well. It should be pointed out that in our TF- 
approximation we did not neglect the centrifugal term l/2r 2 
in Eq. (1241 . which is actually also part of the kinetic energy. 
It is responsible for the vanishing density in the center of the 
vortex core. This turns out to be the better approximation for 
t > 0, while for t < 0, dropping the centrifugal potential 
provides the more accurate approximation because of the ir- 
rotational nature of the initial state. 




FIG. 5: Time evolution of the vortex state density \ip c (r, t)\ 2 vs. ra- 
dius r measured in units of ano- In this simulation the self-scattering 
lengths of 87 Rb were chosen, but at c was deliberately set to zero. 
This change gives rise to an increase of the residual excitations of 
the final vortex state while the transfer efficiency again is higher than 
99%. The dashed and dotted lines show the "breathing" of the vortex 
state at different times t. As a reference, the solid line depicts the 
stationary vortex state. 



B. Total energy of the system 

The total energy of a BE condensed system can be obtained 
from the expectation value of the microscopic Hamiltonian 
with respect to a symmetry -broken ensemble 1371 B8I1 . thereby 
discarding higher order correlation functions. Unless time- 
translational symmetry is broken due to explicit time depen- 
dencies of external fields, the total energy of the system must 
be conserved. Within these assumptions, the energy func- 
tional is given by 



E — E\ao + Ejf + Edi p , 

where the individual contributions are as follows 



(30) 



As shown in Fig. 3, the adiabatic evolution of these sim- 
ple approximations matches the exact numerical results for 



£W = 7r/ drr [\d r ^ a \ 2 + \d r i> b \ 2 + \d r i) c \ 2 } , (31) 
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E TP = iT drr [2A|W| 2 +r- 2 (|^ a | 2 + |V' fc | 2 + |V' c | 2 ) 
Jo 

+ ^\A\ 2 + n(\4, b \ 2 + \^ c \ 2 ) 2 

(32) 

and 

poo 

E Aip = 2ir drr n p ^ + i/>* a rCl d ip c + c.c] . (33) 
Jo 

Eq. ( 13 11 denotes the radial kinetic energy Skin, while Eq. J32i 
represents all the energy arising from the detuning, the trap 
potential, and the mean-field shifts, respectively. The dipole 
energy, which is caused by the dipole coupling of the laser 
fields to the atoms, is given by the expression in Eq. ( 133b . 

The evolution of the total energy and its individual contri- 
butions is shown during the non-equilibrium transfer in Fig.|6] 
As expected, the total energy is constant before and after the 
STIRAP process. However, due to an excitation of the breath- 
ing mode, the kinetic as well as the potential energy exhibit 
complementary oscillations. The dipole energy is negative, 
which is well-known from the interaction of a two-level atom 
with a laser field. This is due to the fact that the polarization of 
the atom is counteracting the external electromagnetic field. 



frequency of these radial excitation is exactly twice the ra- 
dial harmonic trapping frequency, e = 2ui, it must be re- 
lated to the scaling symmetry of the two-dimensional system. 
These "breathing modes" have been studied first theoretically 
in Refs. I43ll44ll and experimentally in Refs. ll45ll . 

The mechanism for exciting these modes arises from 
squeezing the harmonic potential in time. In the present con- 
text of the STIRAP configuration, the origin of such an addi- 
tional potential can be understood from considering Eq. 1241 . 
During the turn-on of the laser fields, we induce an ac-Stark- 
shift potential \fld(r, i)| 2 /A, which is proportional to r 2 . 
While this effect is interesting by itself, it could be eliminated 
easily by an appropriate control of the trapping potential, thus 
reducing the amount of excitations. 

12, . . . n 





FIG. 6: Total energy E (measured in units of Tioj) as a function of 
the time t (in natural units Tko). The different contributions to the 
total energy are shown in the plot. E m =\ denotes the energy of the 
condensate in the stationary vortex state. Because of remaining exci- 
tations the potential and kinetic energy oscillate. However, the total 
energy is conserved before and after the transfer process. 



C. Linear response of the system and breathing modes 

The numerical simulations show that a transfer of almost 
100% of the particles into a vortex state is possible. Still, 
residual radial excitations of the vortex state remain. As the 
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FIG. 7: Mean squared radius (r^)t of the vortex state vs. time (mea- 
sured in units of ano and Tko, respectively). Before the transfer pro- 
cess this state is not occupied. Afterwards the mean squared radius 
oscillates around the equilibrium value eq .) (dashed line). The 
frequency of the excitation is exactly 2uj. 

Moreover, a completely adiabatic transfer is limited by to 
two factors: On the one hand a standard STIRAP process in 
a homogeneous gas requires an adiabaticity condition fir ^ 
1, i.e. the time delay r must be sufficiently large. On the 
other hand r cannot be chose arbitrarily large because here 
we deal with a STIRAP process in a trap. Therefore, coupling 
to external motion has to be taken into consideration for large 
time delays r. 

The residual radial excitations can be visualized by evalu- 
ating the mean squared radius 

2tt r°° 

(r 2 ) t = ^y o drr^ c {r,t)*r 2 ^ c {r,t), (34) 

which is proportional to the potential energy of the vortex- 
component. As shown in Fig. this quantity is oscillating 
with the frequency e = 2 uj, which corresponds to the breath- 
ing mode. 

In addition, small radial excitations (i.e. the linear re- 
sponse) of a BEC can be understood from Bogoliubov the- 
ory lEol l4lll . Therefore, we now consider the radial one- 
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FIG. 8: Mean squared radius (rf)t of the vortex state vs. time (mea- 
sured in units of aHo and Tho, respectively). In this simulation the 
self-scattering lengths of 8 ' Rb were chosen, but at c was deliberately 
set to zero. This change corresponds to an increase of the "breathing" 
of the final vortex state. 




FIG. 9: Excitation frequencies e relative to the ground state fre- 
quency of a quasi-two-dimensional BEC in the vortex state (|m| = 
1) measured in units of u, in comparison to the two-dimensional har- 
monic oscillator with fixed \m\ = 1. For p — and p = 1 the linear 
and nonlinear energies are identical. For p > 2 they split, and the 
excitation energies of the BEC lie below those of the non-interacting 
gas. 



component GP-equation for the vortex-state with angular mo- 
mentum lh that reads 

id t ^(r, t) = (V 1 ' + K |V(r, t)| 2 ) i>(r, t), (35) 

with ht 1 ' defined in Eq. i22\ . With the ansatz 

^(r,t) = e-* 1 (Mr) + u(r)e- l£t + v*(r)e t£t ) , (36) 

where "00 ( r ) is the wave function of the stationary vortex, the 
time evolution of which is determined by the chemical poten- 
tial \i, and u(r) and v(r) denote the excitation modes with the 
normalization (e > 0) 



2tt 



drr (\u(r)\ 2 - \v(r)f 



1. 



(37) 



The spectrum can be calculated from the linear response 
eigenvalue problem 



h 



u(r) 
v(r) 



u{r) 
v{r) 



where 



/ l W- M + 2 K |^ (r)| 2 



(38) 



(39) 



The results of our numerical calculations are shown in 
Fig. [9] and Fig. ^3 being in agreement with earlier work 
14211 . The lowest mode (p — 0) corresponds to the condensate 
wave-function itself (Goldstone-mode). The frequency spec- 
trum can be compared to spectrum of the well-known two- 
dimensional quantum-mechanical harmonic oscillator, which 
is given by 



-HO 



2p- 



1. 



(40) 



where p £ IN denotes the principal quantum number, m E Z 
denotes the angular momentum, and, in our case, |m| = 1 is 
fixed. 
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,y x — ^^r— 


____JD=0 


0.2 
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-0.2 







v(r) 




FIG. 10: Excitation modes obtained by Bogoliubov theory. The first 
three modes of u(r) and v(r) are plotted. The lowest mode (p = 
0) corresponds to the condensate wave-function itself (Goldstone- 
mode). 



IV. CONCLUSIONS 

In this paper, we have developed a novel scheme for the 
optical creation of vortices in a trapped Bose-Einstein con- 
densate, using the technique of stimulated Raman adiabatic 
passage (STIRAP). In our model, we considered a BEC 
of three-level atoms in a A-configuration of the electronic 
states, which are coupled by two co-propagating laser pulses. 
The aim was the transfer of angular momentum, carried by 
one of the beams (Gauss-Laguerre mode), to the BEC. The 
underlying mechanism of STIRAP is analogous to single- 
particle physics. In contrast to the latter case, we derived 
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a multi-component nonlinear Schroedinger equation (Gross- 
Pitaevskii equation), using the mean-field approximation. We 
presented results of numerical simulations that apply to a BEC 
of 87 Rb-atoms. For a suitable set of laser parameters an almost 
100% transfer to the vortex state can be achieved. These re- 
sults can be understood with an intuitive and accurate approx- 
imation within the Thomas-Fermi limit. The occurrence of 
residual radial excitations in the vortex state can be explained 
by so-called breathing modes, which are specific for the two- 
dimensional regime and can be eliminated by an appropriate 



control of the trap frequency. To confirm these explanation, 
we have calculated the Bogoliubov excitation spectrum nu- 
merically. 
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